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Abstract 

We present a simple proof of the factorization of (complex) symmetric matrices 
into a product of a square matrix and its transpose, and discuss its application in 
establishing a uniqueness property of certain antilinear operators. 

1 Introduction 

One of the interesting results of linear algebra is that every square matrix may be factored 
into the product of two symmetric matrices |7j . The factorization of symmetric matrices 
into a product of a square matrix and its transpose is however less known. In fact, there seems 
to be no mention of this factorization in modern texts on linear algebra. The purpose of this 
note is to present a simple derivation of this particular factorization of symmetric matrices 
and to discuss its application in establishing a uniqueness property of certain antilinear 
operators. 
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2 Notation and Definitions 



In this note we shall express the complex-conjugate, the transpose, and the conjugate- 
transpose (adjoint) of a matrix (an operator) m by m, m T , and m*, respectively, and identify 
the elements of C e by columns of i complex numbers. Then the Euclidean inner-product of 
two vectors W\ and W2 takes the form w^Wi- 

Definition 1: Let S be a set, then the Kronecker delta function 5 : S 2 — > {0, 1} is 
defined by 



Definition 2: A function X : 7i — > 7i acting in a complex vector space H. is said to 
be an antilinear operator if for all ijeC and 0, ip G 7i, X{x<f> + y0) = xXcf) + yXip, 
where X(f) means X ((/>). 

Definition 3: An antilinear operator X : H. — > 7i acting in a complex inner-product 
space H with inner-product ( , ) is said to be symmetric or Hermitian if for all 0, ijj 6 Ti, 



Definition 4: A linear operator H : 7i —* Ti acting in a complex inner-product space 
7i is said to have a symmetry generated by a function X : 7i — ► Ti or simply a X- 
symmetry if H and X commute, i.e., [H, X] := HX — XH = 0, where stands for the 
zero operator. A symmetry generated by an antilinear operator is called an antilinear 
symmetry. 

Definition 5: Let H : 7i — > 7i be a linear operator acting in a complex inner-product 
space 7i and G : 7^ — > 7i be an Hermitian invertible linear or antilinear operator. Then 
H is said to be G-Hermitian [|TJ or G-pseudo-Hermitian if H* = GHG^ 1 . 





for a = b 
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3 Motivation: Consequences of antilinear symmetries 

Consider a diagonalizable linear operator H : 7i — > 7i acting in a finite-dimensional complex 
inner-product space 7i with inner-product ( , ). Let n label the eigenvalues E n of H, fi n 
be the multiplicity of E n , and i[) n)a be the eigenvectors corresponding to the eigenvalue E n 
where a G {1, 2, • • • , //„} is the degeneracy label. Then it is well-known || that the adjoint 
H* of if is diagonalizable; the eigenvalues E n of H* are complex conjugate of those of H, i.e., 
E n — E n ; the multiplicity of _E n is equal to fi n ; and one can choose the eigenvectors 0„ ia of H* 
in such a way that for all spectral labels m, n and degeneracy labels a, b, (<p n ,a, ipn,a) — 8n,m$a,b- 
Clearly, both sets of eigenvectors ipn,a of H and eigenvectors ai „ of H* form bases of H; 
{ip n ,ai <Pn,a} is a complete biorthonormal system. Recently 0, we have shown that if the 
eigenvalues of H are real, then H has an antilinear symmetry. More generally, we proved 
the following theorem. 

Theorem 1: The presence of an antilinear symmetry of H is a necessary and sufficient 
condition for the eigenvalues of H to either be real or come in complex conjugate pairs. 

The proof of Theorem 1 uses the following lemma. 

Lemma 1: Every diagonalizable linear operator H : TC — ► 7i acting in a finite- 
dimensional complex inner-product space Ti, is T-Hermitian, 

H* = THT-\ (1) 

for some Hermitian, invertible, antilinear operator T : Ti — ► TC. 

It turns out that any such T may be expressed in terms of the eigenvectors 0„ ia of H* 
according to 

V(eH, T( = E °ba (^a. O0n,6, (2) 

n a,b=l 

where are the entries of symmetric invertible /i n x /i n matrices . 
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Note that Theorem 1 and Lemma 1 have infinite-dimensional generalizations for linear 
operators H admitting a complete biorthonormal system of eigenvectors ||. 
Next, consider a general basis transformation, 



where are the entries of an invertible fi n x fi n matrix v^ n \ In terms of the transformed 
basis vectors T has the form: V£ G TC, T( = n b=1 c^\<f/ na , C)4>' n &, where c'^ are 

related to according to = v^c'^v^ 7 '. This equation indicates that the issue of the 
uniqueness of T for a given H is related to whether one can find for each an invertible 
matrix such that = v^ n 'v^ n > T . In the remainder of this note we shall give a proof of 
the fact that this is indeed possible, and one can transform to a basis in which T has the 
(canonical) form (fj) with c^j = 5 a b for all n. 

4 Factorization of Symmetric Matrices 

Theorem 2: A square matrix c is symmetric if and only if it can be written as c = vv T 
for some square matrix v. 

Proof: If c = vv T then clearly c is symmetric. To prove the converse we use induction 
on the dimension n of the matrix c. For n = 1, c = vv T = v 2 is trivially satisfied by 
letting v := yfc. By induction hypothesis we assume that for all k e {2, • • • , n}, every 
k x k symmetric matrix c can be written in the form c = vv T for some k x k matrix 
v. Now let C be an (n + 1) x (n + 1) symmetric matrix. Then C has at least one 
eigenvector M i.e., there are e G C n+1 — {0} and A G C such that 



Now let V := {w G C n+1 |w*e = 0} be the orthogonal complement of e. Clearly V is an 
n-dimensional vector subspace of C n+1 . Next, consider the following two possibilities. 




(3) 



6=1 



Ce = Ae. 



(4) 
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(i) e ^ V. In this case, choose a basis {ei, e n , ■ • • , e„} of V and let e n+ i := e. Then 
e*2, • • • , e n , e n+ \} is a basis of C n+1 and the matrix A := (ei, e n , ■ ■ • , e n , e n+ i) 
is invertible. Note that for all I G {1, 2, • • • , n}, 6 V, and efe n+ \ = 0. This in 
turn implies that the matrix A T CA which is symmetric has the block form 



A T CA 



c 
T A 2 



(5) 



where c is a symmetric n x n matrix. By induction hypothesis there is an n x n 
matrix v such that c = vv T . Now let B be the (n + 1) x (n + 1) matrix 

S T 



T A 



(6) 



and V := (5A- 1 ) T . Then in view of (g) and (§), = A T CA and 



This completes the proof for case (i). 

^ii) e G V, i.e., e T e = 0. In this case, let V := {to G V|?Z>*e = 0} be the or- 
thogonal complement of e in V, {ei, e*2, • ■ ■ , e n _i} be a basis of V, e n = e, and 
e„ + i = e. Then {ei,e 2 , ••• , e n ,e„ + i} is a basis of C" +1 and the matrix A' : = 
(e*i, e n , • • • , e n , e n+ i) is invertible. Note that for all t G {1, 2, • • • , n — 1}, e^ G V, 



and e/"e n +i = 0. Furthermore, e^ +1 e n+ i = e T e = and a := eje n+ i = e*e G 



In view of these relations and 

/ 



C := A' T CA' 



~<1 


C l,2 


C 'l,n-1 


C l,n 


^2,1 


^2,2 


^2,n-l 


^2,n 



n-1,1 Wi-1,2 



V 



C n,l 




C n,2 







n— l.n— 1 n— l,n 



r' 







Aa 



\ 





Aa 
j 



(7) 



where c[ ^ := e?Cej are the entries of a symmetric nx n matrix c'. Now if A = 0, 
C is block-diagonal and the argument given in case (i) leads to a proof of the 



theorem. This leaves the case A 7^ 0. In this case, let A := A'D where D is the 
(n + 1) x (n + 1) matrix 

( 1 ••• x x \ 

1 ••• x 2 

D 

-•• 1 i n _ :l 

■■■ 1 x r 

\Vl Vl VZ ■■■ Vn-l y n 
with Xi, X2, ■ • ■ x n -i being arbitrary complex numbers, 



Xv 



(9) 



and for all % e {1, 2, • • • , n} 



(10) 



Then a simple computation shows that 



A T CA = D T C'D = [ C J J ° 
\ T A 

where the entries of a symmetric nxn matrix c and A' G C. Therefore, we 

can use the argument given in case (i) to show the existence of an (n+1) x (n+ 1) 
matrix B satisfying 



AC A = B B. 



The proof of the theorem will be complete if we show that the matrix A = 
A'D is invertible. Because A' is invertible, it suffices to show the existence of 
X\, X2, ■ ■ ■ x n -\ for which det D^O. We can use the properties of the determinant 
and Equations @ and fllCf ) to compute 



det D = - XiUi 
i=i 

n—1 n—1 n—1 

i=l i<j=l i=l 



Suppose that for all values of x±, x 2 , 

in which case C will have the form 

/ 




c 



■x n _i, det D = 0. This implies that d = 



•0 \ 
•000 



■•■ 
■■• Act 
••• Xa 



By induction hypothesis we have a 2 x 2 matrix m satisfying 




m T m. 



Therefore, setting 



/ 




B :-- 














\ 





• 

• 777,1,2 ?77i >2 

V • • • m 2 ,i m 2j 2 J 

and V := {BX- l ) T , we have C = B T B and C = VV T . If c' ^ then there 
are values of xi,x 2 , ■ • ■ x n -\ for which D is a nonsingular matrix and A = A'D 
is invertible. Hence we can set V := (BA- 1 ) 7 and use QTT]) to show that C = 
VV T . □ 



Concluding remarks 

1. The factorization c = vv T established in Theorem 2 is invariant under the transforma- 
tion v — > f ' = vo where o is an arbitrary (complex) orthogonal matrix. In particular, 
one may choose o so that the factorizing matrix v' has a simple form. 

2. In view of the discussion of Section 3, one has the following consequence of Theorem 2. 



Corollary 1: Up to basis transformations ([|), there is a unique antilinear oper- 
ator T satisfying Equation (fj), namely 

/-in 

VC G H, n = J2 E^> O0n,a, (13) 

n a=l 

3. For a self-adjoint linear operator H, one can set 0„ iQ = ip n>a and use the completeness 
of the eigenvectors ip n a and (|T3"D to deduce T 2 = /, where / is the identity operator. 



Furthermore, noting that in this case ([I]) is equivalent to T-symmetry of H, one can 
prove the following. 

Corollary 2: Every self-adjoint linear operator H has an antilinear symmetry 
generated by a Hermitian, invertible, antilinear operator T satisfying T 2 = I. 
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